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Abstract. Some BL Lacertae objects show a periodic 
behaviour in their light curves that is often attributed 
to the orbital motion of a central binary black hole sys- 
tem. On this basis, and assuming a circular orbit, Rieger 



and Mannheim (2000) have recently proposed a method 



to determine the orbital parameters of the binary system 
from the observed quantities, i.e. the signal periodicity, 
the flux ratio between maximum and minimum signal and 
the power law spectral index of the photon flux. 

However, since these binary black holes are expected 
to originate from galactic mergers, they could well be on 
eccentric orbits, which might not circularize for a substan- 
tial time. We therefore gen eralize the treatment proposed 
by Rieger and Mannheim ( 2000| ) by taking into account 
the effect of the orbital eccentricity of the binary system. 
We apply the model to three well-observed Markarian ob- 
jects: MKN 501, MKN 421 and MKN 766 that most likely 
host a binary system in their centers. Some astrophysical 
implications of this model are also investigated with par- 
ticular emphasis to the gravitational radiation emission 
from the binary black holes. Under particular conditions 
(e.g. for some values of the orbital separation, eccentric- 
ity and Lorentz factor) one can obtain signals above the 
sensitivity threshold of the LISA detector. 
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1. Introduction 

BL Lacertae objects, also known as Markarian objects 
(hereafter MKNs), belong to the class of active galaxies. 
According to the well-established unified model on radio- 
loud active galactic nuclei (see Urry and Padovani 1995| ), 
these objects are thought to be dominated by relativistic 
jets seen at small angles to the line of sight. The structure 
of the relativistic jets is largely unknown and the typical 
smooth and fairly featureless spectra of such objects may 
be reproduced by models making use of very different as- 



sumptions. Combining spectral and temporal informations 
greatly constrains the available models. 

Up to now, several astrophysical phenomena have been 
attributed to bi nary b lack holes, like precession of jets 
(Be gelma n et al. 1980 ), misalignment (Conway and Wro- 
bel |1995| ), periodic outburst activity in the quasar OJ 288 
(Sillinpaa et al. 1988, Letho and Valtonen |1996 ) and pre- 
cession of the accretion disk under gravitational torque 
(Katz |l997p . 

It has been recently observed that some MKN objects 
show a periodic behaviour in the radio, optical, X-ray and 
7-ray light curves that is possibly related to the presence 
of a massive binary black hole creating a jet either aligned 
along the line of sight or interacting with an accretion disk 
(Yu |200l|) . Therefore, the search for X-ray /7-ray variabil- 
ity can be considered as a method to probe the existence 
of a massive binary black hole in the center of a galaxy. 

Massive black hole systems are expected to be fairly 
common in the Universe as a result of merging between 
galaxies. Indeed, in the framework of the unified model 
for th e morphological evolution of galaxies (see e.g. White 
1997 ) giant ellipticals, which are the host galaxies of most 
MKN objects, appear to be the result of mergers between 
spirals. On the other hand, it is now well established that 
galaxies g enera lly contain massive black holes i n their nu- 
clei ( Rees 1984 , Kormendy and Richstone 1995 , Richstone 
et al. 1998 ). Therefore, merging would naturally lead to 



the form ation of massive binary black holes (Begelman et 
al. [19801 ). 

At least three MKN objects (i.e MKN 501, MKN 421 
and MKN 766) are particularly well studied at high ener- 
gies, revealing a possible periodic behaviour in their light 
curves. 

MKN 501, at z = 0.034, shows a clear well-correlated 
23 day periodicity in X-ray and TeV energy bands with an 
observed TeV flux ratio / — 8 between the maximum and 
minimum of the signal (Protheroe et al. 1998] , Hayashida 



et al. |1998j , Kranich et al. |1999[ Nishikawa et al. |1999| ), 



while evidence for correlations in the optical U-band is 



rather weak (Catanese et al. 1997 , Aharonian et al. 1999) 
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It has also been suggested that the complex morphology 
of the jet and the peculiar behaviour of its spectral energy 



2 



F. De Paolis et al.: Astrophysical implications of binary black holes 



distribution are probably related to the presence of a mas- 



sive binary black hole (Conway and Wrobel 1995, Villata 
and Raitcri |1999| ). 

MKN 421, at z = 0.031, is the brightest BL Lacertae 
object at X-ray and UV wavelengths and it is the first 
extragalactic source discovered at TeV energies (Punch et 
al. 1992). This nearby source, which has been recently ob- 
served by the XMM-Newton (see Brinkmann et al. [2000 ) 



and by BeppoSax (Maraschi et al. 1999) satellites, shows 
remarkable X-ray variability correlated with strong activ- 



ity at TeV energies (George et al. 198S) on a time-scale of 
~ 10 4 s (Maraschi et al. |1999 ) and with a flux ratio / — 2. 

X-ray observations of the nearby MKN 766, at z = 
0.013, have been performed by the XMM-Newton satel- 
lite (Boiler, Keill, Tritmper et al. 2000). These observa- 
tions have revealed the presence of a strong X-ray peri- 
odic signal with frequency ~ 2.4 x 10 -4 Hz and flux ratio 
/- 1-3. 

Based on the assumption that the periodic observed 
light curve of MKN 501 is related to the presence of a 
binary system of black holes (one of which emits a jet), 
Rieger and Mannheim ( 2000| ) have proposed a method to 
determine the physical parameters of the binary system 
from the observed quantities, i. e. signal periodicity, flux 
ratio between maximum and minimum signal and power 
law spectral index. They adopted the simplifying assump- 
tions that i) the two massive black holes move on circular 
orbits around the common center of mass and ii) the bi- 
nary separation a is such that the gas dynamical time-scale 
T gas is equal to the time-scale for gravitational radiation 



T gw (see also Begelman, Blandford and Rees 1980) 



However, if binary black holes originated from galactic 
mergers, they could be on eccentric orbits and eccentricity 
values up to 0.8-0.9 are not necessarily too extreme (Fitch- 
ett 1987). Of course, due to gravitational wave emission, 
orbits tend to circularize but this happens within a time- 
scale of the same order of magnitude as the merging time- 



scale (Peters 1964, Fitchett 1987). Therefore, if a massive 



binary black hole is found at the center of a galaxy, it 
may happen that the constituting black holes are still on 
eccentric orbits, in whi ch ca se the method proposed by 
Rieger and Mannheim ( 200C ) does not hold. 

As far as point ii) is concerned, we note that the as- 
sumption that nowadays T gas is equal to T gw is arbitrary. 
Indeed, it is possible that T gas = T gw at some time in 
the past, but since then the binary system evolution was 
driven by gravitational wave emission which shrunk the 
orbital separation faster (see Section 2). 

The aim of the present paper is to generalize the model 



proposed by Rieger and Mannheim (2000) by taking into 
account the effects of the orbital eccentricity of the binary 
system and considering the binary separation a as a free 
model parameter. This is done in Section 2, where the 
model is also applied to three well-observed MKN objects 
(MKN 501, MKN 421 and MKN 766) obtaining a set of 
possible values for the free model parameters jb, o, and e 



corresponding to different values for the masses m and M 
of the two black holes. In Section 3 we study the gravita- 
tional waves emitted by the binary systems at the center 
of the above-quoted MKN objects and analyze the possi- 
bility of detecting them by the LISA gravitational wave 
interferometer. Our conclusions are presented in Section 
4. 

2. Binary black hole model for BL Lacertae 
objects 

In the commonly-accepted evolutionary scenario, giant el- 
liptical galaxies are believed to be the product of mergers 
between spiral galaxies. Since each galaxy likely contains 
a central black hole with mass 10 6 -10 9 M (Richstone, 



Ajhar and Bender [1998 ), merging would sometimes pro- 
duce binary black hole systems separated by a typical dis- 
tance of 0.1 — 1 pc (Begelman, Blandford and Rees 1980| )). 

We assume that binary black holes exist in the center 
of many MKN objects and in particular in those of MKN 
501, M KN 421 and MKN 766. As in Rieger and Mannheim 
( 2000| ) , the periodicity in the flaring state, observed to- 
wards these objects, is assumed to be the consequence of 
the orbital motion of a relativistic jet in the binary black 
hole. Therefore, the observed signal periodicity has a ge- 
ometrical origin, being a consequence of Doppler-shiftcd 
modulation. Let us further assume that the jet is emitted 
by the less massive black hole and that the nonthermal 
X — ray /7-ray radiation propagates outwards from the 
core along the jet with Lorentz gamma factor 7;,. The ob- 
served flux modulation due to Doppler boosting can be 
written as 



S(v) = 5 d+a S'(v) , 



(1) 



where a is the source spectral index Q and the Doppler 
factor is given by 



5 = 



Vl-K 2 + ^ 2 s )/c 2 
1 — (v z cos i + vis sin i)/c 



(2) 



Here, v z is the outflow velocity in the direction of the total 
angular momentum, i is the inclination angle between the 
jet axis and the line of sight and vi s is the component of 
the less massive black hole velocity along the line of sight. 

Let M and m represent the masses of the primary and 
the secondary black hole, respectively. Conventionally, the 
problem of predicting the motion of a system of two grav- 
itationally interacting point masses is simplified by con- 
sidering the equivalent system. This consists of a mass 
M t = M + m fixed in space and acting on a reduced mass 
/i = Mm I (M + m) that orbits around it. We assume that 



1 Values of the power law index a for the three MKNs of 
interest are found to be 1.2, 1.7 and 2.11 for MKN 501, MKN 
421 and M KN 766 , respectively. For more details see Riege r and 
Mannheim (2000), Guainazzi , Vac anti, Maizia et al. ( 199S| ) and 
Boiler, keill, Triimper et al. (200C|). 
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the reduced mass moves around the total mass on an el- 
liptic orbit with semi-major axis a and eccentricity e. Let 
(r, 9) be the polar coordinates of (x relative to M t , so that 
the equation of the orbit is given by 



a(l-e 2 ) 
1 + e cos 9 



(3) 



It is straightforward to demonstrate that the component 
of the reduced mass velocity along the line of sight is given 
by (Smart $97% 



IG{M + m) 
o(l - e 2 ) 



sin f 



(4) 



The component of the less massive black hole velocity 
along the line of sight is thus given by 



M 



in 



+ m 

which can be rewritten as 



Vis 



Ra^-a 



sin# 



(1-e 2 )! 



(5) 



(G) 



where R a = Ma/(M + m) and Vl a = [G(M + m)/a 3 ] 1/2 
is the Keplerian orbital frequency if the mass m moves 
on a circular orbit of radius a around M. The velocity 
vis is maximal for 9 = ir/2 and minimal for 9 — 3tt/2, 
corresponding through equation (^) to the two possible 
values S max and 5 m i n of the Doppler factor. 

With the assumption that the periodicity in the ob- 
served signal is due to the orbital motion of the binary 
black hole, from equation (Q) one obtains the condition 
Smax/S m in — f 1 ^ 3+a \ where / is the observed maximum 
to minimum flux ratio. Consequently, by using equations 
(Q) and (||) we have 



Q a Ra 



f 



1/3+Q 



yl/3+a _|_ 1 y s [ n j 



cotz (1 - e ) 2 c 



(7) 



which, in the limit e — > 0, reduces to equation (4) in Rieger 
and Mannheim (|2000D . 

As a consequence of the binary orbital revolution 
around the center of mass, the observed signal period P a b s 
is related to the keple rian period P^ — 27r/O a by (Rieger 
and Mannheim, 2000) 



(8) 



P obs = {± + z) (l-^cosz)p fc , 

where i ~ l/j b and v z ~ c(l - 1/7 2 ) 1/2 (Spada |199' . 
From equations (0) and (||), one obtains the following re- 
lation 



M 



P. 



1/3 



(m + M) 2 / 3 [2tt(1 + z)G] 1 / 3 



j-l/3+a _ j 
JTJ3+ 



(9) 



, 1 - 
1 V c 



— COS I J (1 

c / 



which, in the limit 



Rieger and Mannheim (2000) 
Moreover, from equation ( 



0, reduces to equation (8) of 
one has 



M = 



2tt(1 + z) (1 



COS I 



-i 2 



Po 



hs 



a 

~G 



(10) 



where a represents the semi-major axis of the binary sys- 
tem. 

Therefore, the allowed primary and secondary black 
hole masses m and M are obtained by solving simultane- 
ously equations (||) and (0). In these equations 7b, e and 
a have to be considered as free model parameters to be 
determined by the observed data (P b s , f and a) towards 
the three MKK of interest. In this way, we can estimate 
the masses of the black holes supposed to be at the centers 
of MKN 501, MKN 421 and MKN 766. 

We mention that equation ( |l0| ) is different from the 
corresponding one in Rieger and Mannheim ( [2000 ) , where 
the binary separation a is obtained by equating the gas dy- 
namical time-scale (Begelman, Blandford and Rees 1980) 
with the time-scale for gravitational radiation (Peters and 
Mathews 1963). Here, as previously stated, we release this 



assumption and consider a as a free fit parameter. 

For the three MKN objects considered, the secondary 
mass m as a function of the primary one M is shown in 
Fig. |] for selected values of -fb = 10, 15 and e = 0.5. For 
each binary system, the intersection between lines with 
the same Lorentz factor determines the masses of the two 
black hole components. 

In Fig. (for MKN 501) we show the primary and 
secondary black hole masses as a function of the orbit 
eccentricity for two values of the Lorentz factor jb = 10, 30 
and a = 8 x 10 16 cm. In this case, the black hole masses 
are in the range 10 6 — 10 9 Mq. 



3. Gravitational waves from massive black hole 
binaries in Markarian objects and comparison 
with LISA sensitivity 

In this Section, the massive black hole binaries in the 
MKN objects considered in Section 2 are examined as pos- 
sible sources of gravitational radiation. Due to the orbital 
period of the binary systems considered, it is expected 
that they emit most of the gravitational radiation in the 
low- frequency band {y < 10~ 4 Hz). Gravitational waves in 
this band are only accessible to experiments based on the 
Doppler tracking technique of inter-planetary spacecrafts. 
The Laser Interferometer Space Antenna (LISA), consist- 
ing of a constellation of three drag-free spacecrafts at the 
vertices of an ideal equilater triangle with sides of ~ 5 x 10 6 



km (Hiscock, Larson, Routzahn and Kulick 2000), has an 



optimal sensitivity in the frequency rang e covering the 
band 10~ 5 Hz - 3 x 10" 2 Hz (Vecchio |l999|) . LISA will be 
able to carry out a deep and extensive census of black hole 
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Fig. 1. From the upper to the bottom panel, the required 
mass dependence for binary black hole model is shown in 
the case of MKN 501, MKN 421 and MKN 766, respec- 
tively. The thick and thin lines represent the conditions 
expressed in equations (§|) and (|l^) with Lorentz factor 
7b = 10 and jf, = 15. The binary separation a is set to 
5 x 10 16 cm, 2 x 10 14 cm and 1 x 10 13 cm for MKN501, 
MKN421 and MKN766, respectively. We have set the orbit 
eccentricity value equal to 0.5. The intersection between 
lines corresponding to the same Lorentz factor gives the 
masses of the black holes in the binary system considered. 



populations in the Universe, providing an accurate demog- 
raphy of these objects and their interactions. In particular 
we expect that binaries of massive black holes could be de- 
tectable by LISA even if the emitted gravitational waves 
are characterized by extremely low frequencies (y < 10 
Hz). This could be the case for the binaries at the cen- 
ter of the three MKN objects considered under particular 
conditions (e.g. for the same values of a, e and 7b). 

Indeed we expect that in a merging process between 
galaxies, the two colliding black holes at their centers 
would eventually move on orbits with non-negligible ec- 
centricity. This fact should have great importance as far 
as the detection of the emitted gravitational waves is con- 



Fig. 2. The dependence of the primary and secondary 
mass on the eccentricity of the orbit is shown for MKN 
501 Lorentz 7^ factors 7& = 10 (dotted lines) and 7t, = 30 
(solid lines), respectively. Here we have set a — 8 x 10 16 
cm. Similar plots can be obtained for the MKN421 and 
MKN766. 



cerned. In fact, a binary system with an eccentric orbit 
emits significantly more power in gravitational waves than 
a circular system with the same semimajor axis. Moreover, 
eccentric binaries will also emit gravitational waves in a 
multitude of harmonics of the orbital frequency u>k associ- 
ated with the semimajor axis a. Thus, an eccentric system 
emits gravitational waves at the frequencies u) n = nu>k 
(with n — 1,2,3...), whereas circular binaries emit purely 
in the n = 2 mode. Peters and Mathews ( 1963| ), found that 
the orbit averaged power emitted in gravitational waves 
by a binary system with eccentricity e in the nth harmonic 
of the orbital frequency is 



dE(n,e) _ 32 G 4 M 2 m 2 (M +■ 



dt 



-g{n,e) , 



(11) 
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where the function g(n, e) depends on n and e by the 
following relation 

4 

71 

g(n, e) = — {[J n ^ 2 {ne] - 2eJ n _i(ne) 



H — J nine) + 2eJ n+ i(ne) - J n+2 (ne) 



+ (1 - e 2 ) [J n _ 2 (ne) - 2J„(ne) + J„ +2 (ne)]' 



(12) 



3n 2 



[J„(ne)]' 



Here J„ (with n integer) is the usual Bessel function of 
order n th . The characteristic spectral amplitude at the 
frequency ui n can be evaluated as 



h{u) n ) 



TTUlnD 



G dE(n, e) 



dt 



(13) 



where D is the distance of the gravitational wave source 
from Earth. 

Equations (|Il|)-(|l3|) allow us to evaluate the expected 
spectrum of the gravitational waves emitted by the mas- 
sive black hole binaries assumed to be in the center of the 
considered MKN objects. Assuming 75 = 5, Fig. || shows 
the expected gravitational wave amplitude h as a function 
of the frequency w, for different values of the binary orbit 
eccentricity. In each panel of the Figure we give for com- 
parison the LISA sensitivity curve (the oblique solid line) 
for an integration time of 5 years. As can be noted, in sev- 
eral cases the expected gravitational waves emitted by the 
considered binary systems may be detected by LISA satel- 
lite. Eccentricity values up to 0.8 are not necessarily too 
extreme. In fact, if a black hole binary forms in a merging 
process, a very high value of eccentricity is expected. Of 
course, due to the gravitational wave emission the orbit 
tends to circularize but this should happen within a time- 
scale of the same order of the magnitude of the merging 
time-scale (Peters ^964). 

In Fig. ^, the binary semimajor axis a is shown as 
a function of the orbital eccentricity e and for two dif- 
ferent Lorentz factors j b for MKN 501, MKN 421 and 
MKN 766, respectively. Each point in the 3D scatter 
plot fits the observations and satisfies eqs. @ and (|l(J). 
Each point on the plot corresponds to a different distance 
d = h(u) n )/hLisA(k>n) (vertical axis) between the gravi- 
tational wave emitted spectrum and the LISA sensitivity 
curve for an integration time of 5 years (the minimum 
distance between the curves diminishes as d increases). 
In Fig. ^, only models with infd(u; n ) > 0.1 are shown. 
Models corresponding to d > 1 in the plot have to be con- 
sidered as observable by the LISA satellite. We notice that 
models corresponding to higher values of the parameter d 
have typically higher masses and lower coalescing times 
with respect to those with lower d values. In any case, 



4. Discussion 

We have considered some astrophysical implications of the 
existence of massive black hole binary systems at the cen- 
ter of MKN objects. These systems are expected to be 
fairly common in the Universe as a consequence of merg- 
ing processes between galaxies. In fact, following the evo- 
lutionary scenario, giant elliptical galaxies are thought to 
be the result of merging between spirals. Since almost each 
galaxy contains a massive black hole in its center, merging 
process would naturally lead to the formation of massive 
black hole binaries. 

Currently, the only method to probe the existence of 
massive binary black holes is to search for double nuclei 
in the galactic c enters . Howev er, as pointed out by Rieger 
and Mannheim ( 2000 ) and Yu ( 2001 ), a periodic behaviour 
in the observed radio, optical, X-ray or 7-ray light curves 
could be considered as the signature of a massive binary 
black hole. 

In this paper, we focused on the possibility of explain- 
ing the intensity variations observed in some MKN objects 
as the consequence of the motion of binary black holes 
around the common center of mass. We have generalized 
the treatment of Rieger and Mannheim ( 200C ) considering 
in detail the effects of the orbital eccentricity. 

At least three MKN objects (MKN 501, MKN 421 and 
MKN 766) reveal a possible periodic behaviour in the ob- 
served signal. In particular, it has been recently addressed 
whether MKN 501 hosts a massiv e black hole binary sys- 
tem (Rieger and Mannheim 2000] ). The same likely holds 
for MKN 421 and MKN 766. The model described in Sec- 
tion 2 allows us to estimate the mass of each black hole in 
the binary system (see Fig. |l|) as a function of the Lorentz 
factor 7b, of the orbital radius a and eccentricity e. 

If massive binary black holes really exist in the cen- 
ters of the MKN objects considered above, they should 
also emit low frequency gravitational waves typically with 
lu < 10~ 4 Hz. In Fig. |] we have presented the gravitational 
wave spectrum emitted by binaries that fit all the obser- 
vational data (flux ratio /, the signal periodicity P b s and 
the source spectral index a). As one can note, for some val- 
ues of a and e one can obtain curves above the sensitivity 
threshold of the LISA detector (we assume an integration 
time of 5 years). 

It is obvious that binary systems emitting gravitational 
waves detectable by the LISA satellite have short coalesc- 
ing times. For example, in the case of MKN 501 the av- 
erage coalescing time for the models represented in Fig. 
[I| (for 7f, = 5) is ~ 10 3 years and tends to decrease for 
higher values of the parameter d. 

However, for signal detection, some caution is needed 
since the LISA sensitivity curve we use in Figs. |^ and 
I] does not take into account the noise from the galactic 
binaries containing two compact objects. Especially im- 
portant are galactic white dwarf binaries (Hils and Ben 
dcr 2000 Nelemans, Yungelson and Portegies-Zwart 2001 
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Fig. 3. The amplitude of the gravitational waves emitted by the binary black holes fitting all the observational data is 
reported as a function of the frequency for the three MKN objects considered. The continuous lines from left to right 
in each panel correspond to orbital eccentricity e = 0.5, 0.6, 0.7, 0.8, respectively (for MKN 501 we added a curve 
corresponding to e = 0.85). The chosen orbital radius a is indicated in each panel. The solid oblique line represents 
the LISA sensitivity with an integration time of 5 years. 



Nelemans, Portegies-Zwart and Verbunt 1999Q , primordial 
black holes (Hiscock 1998) and white dwarf binaries (His- 
cock, Larson, Routzahn and Kulick 2000) in the galactic 



halo. These binaries would give rise to a confusion noise 
which lies well above LISA's sensitivity curve in the fre- 
quency range of interest between 10 -4 Hz and 10~ 3 Hz 
(see e. g. Fi g. 3 in Nelemans, Yungelson and Portegies- 



Zwart 2001 ). This noise would make even more difficult 
the detection of the considered massive black hole binary 
systems so that, in order to have detectable signals, one 
would need more extreme conditions (e.g. higher values of 
the parameter d or, equivalently, shorter coalescing times). 



As a final remark we notice that it is possible that 
there are a lot of hard binary black holes at the center of 
galaxies, MKN objects, Seyfert galaxies, etc., and that we 
can observe directly (from the periodicity of lightcurves at 
high energy) only a small fraction of them. If this is so, de- 
tection of gravitational waves from such systems would be 
a means to detect them. It is worth noting that the num- 
ber of intervening galaxies between Earth and a galaxy at 
redshift z is 

2 



N ~ 6x10" 



R 



N„ 



r ) [{l+zf' 2 -l] ,(14) 



v 3kpcy VO.lMpc" 

where R is the typical galaxy radius and N g the average 
number density of galaxies. Therefore, it is expected that 
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Fig. 4. In the 3D scatter plot, the binary semimajor axis a is shown as a function of the orbital eccentricity e for 
different values of Lorentz factor 75, for MKN 501, MKN 421 and MKN 766, respectively. Each point in the evidenced 
regions of the a-e plane fits the observations and satisfies eqs. (^|) and (|o|). On the vertical axis, the detection parameter 
d is also shown (for detalis see text). 



at least in some cases gravitational lenses, displaced along 
the line of sight of a gravitational wave source, should dis- 
tort the path of the gravitational waves just as they dis- 
tort the path of the observed light. For gravitational waves 
with frequency v < 10 -4 Hz, the LISA angular resolution 



is below 10~ 2 steradians (Cutler 1998 ). Consequently, the 
lensed source is always seen as an amplitude amplification 
(i.e. as a microlensing event) and the observed amplitude, 
at the frequency cu n , is (De Paolis, Ingrosso and Nucita 
2001]) 



h max (uj n ) ~ VA=/iK) 



(15) 



where A max is the amplification parameter depending on 
the geometry of the microlensing configuration. This am- 
plification could be in some cases of the order of 10 or 
more, thereby increasing the number of detectable binary 
black holes at galactic centers by future gravitational wave 
detectors like LISA. Even more interesting would be the 
effect of the gravitational wave time-delay due to the dif- 
ferent paths traveled by the deflected signals. This time- 
delay, of about one year as for the lensed Quasars, would 
give rise to a frequency difference due to the binary evolu- 
tion with time. This effect could produce a "beating pat- 
tern" in the wave amplitude which could be searched for 
in future gravitational wave experiments. 
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Abstract. Some BL Lacertae objects show a periodic 
behaviour in their light curves that is often attributed 
to the orbital motion of a central binary black hole sys- 
tem. On this basis, and assuming a circular orbit, Rieger 
and Mannheim (2000) have recently proposed a method 
to determine the orbital parameters of the binary system 
from the observed quantities, i.e. the signal periodicity, 
the flux ratio between maximum and minimum signal and 
the power law spectral index of the photon flux. 

However, since these binary black holes are expected 
to originate from galactic mergers, they could well be on 
eccentric orbits, which might not circularize for a substan- 
tial time. We therefore generalize the treatment proposed 
by Rieger and Mannheim (2000) by taking into account 
the effect of the orbital eccentricity of the binary system. 
We apply the model to three well-observed Markarian ob- 
jects: MKN 501, MKN 421 and MKN 766 that most likely 
host a binary system in their centers. Some astrophysical 
implications of this model are also investigated with par- 
ticular emphasis to the gravitational radiation emission 
from the binary black holes. Under particular conditions 
(e.g. for some values of the orbital separation, eccentric- 
ity and Lorentz factor) one can obtain signals above the 
sensitivity threshold of the LISA detector. 

Key words: BL Lacertae objects: general - Black hole 
physics 



1. Introduction 

BL Lacertae objects, also known as Markarian objects 
(hereafter MKNs), belong to the class of active galaxies. 
According to the well-established unified model on radio- 
loud active galactic nuclei (see Urry and Padovani 1995), 
these objects are thought to be dominated by relativistic 
jets seen at small angles to the line of sight. The structure 
of the relativistic jets is largely unknown and the typical 
smooth and fairly featureless spectra of such objects may 
be reproduced by models making use of very different as- 

Send offprint requests to: F. De Paolis 



sumptions. Combining spectral and temporal informations 
greatly constrains the available models. 

Up to now, several astrophysical phenomena have been 
attributed to binary black holes, like precession of jets 
(Begelman et al. 1980), misalignment (Conway and Wro- 
bel 1995), periodic outburst activity in the quasar OJ 288 
(Sillinpaa et al. 1988, Letho and Valtonen 1996) and pre- 
cession of the accretion disk under gravitational torque 
(Katz 1997). 

It has been recently observed that some MKN objects 
show a periodic behaviour in the radio, optical, X-ray and 
7-ray light curves that is possibly related to the presence 
of a massive binary black hole creating a jet either aligned 
along the line of sight or interacting with an accretion disk 
(Yu 2001). Therefore, the search for X-ray/7-ray variabil- 
ity can be considered as a method to probe the existence 
of a massive binary black hole in the center of a galaxy. 

Massive black hole systems are expected to be fairly 
common in the Universe as a result of merging between 
galaxies. Indeed, in the framework of the unified model 
for the morphological evolution of galaxies (see e.g. White 
1997) giant ellipticals, which are the host galaxies of most 
MKN objects, appear to be the result of mergers between 
spirals. On the other hand, it is now well established that 
galaxies generally contain massive black holes in their nu- 
clei (Rees 1984, Kormendy and Richstone 1995, Richstone 
et al. 1998). Therefore, merging would naturally lead to 
the formation of massive binary black holes (Begelman et 
al. 1980). 

At least three MKN objects (i.e MKN 501, MKN 421 
and MKN 766) are particularly well studied at high ener- 
gies, revealing a possible periodic behaviour in their light 
curves. 

MKN 501, at z = 0.034, shows a clear well-correlated 
23 day periodicity in X-ray and TeV energy bands with an 
observed TeV flux ratio / — 8 between the maximum and 
minimum of the signal (Protheroe et al. 1998, Hayashida 
et al. 1998, Kranich et al. 1999, Nishikawa et al. 1999), 
while evidence for correlations in the optical U-band is 
rather weak (Catanese et al. 1997 , Aharonian et al. 1999). 
It has also been suggested that the complex morphology 
of the jet and the peculiar behaviour of its spectral energy 
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distribution are probably related to the presence of a mas- 
sive binary black hole (Conway and Wrobel 1995, Villata 
and Raiteri 1999). 

MKN 421, at z = 0.031, is the brightest BL Lacertae 
object at X-ray and UV wavelengths and it is the first 
extragalactic source discovered at TeV energies (Punch et 
al. 1992). This nearby source, which has been recently ob- 
served by the XMM-Newton (see Brinkmann et al. 2000) 
and by BeppoSax (Maraschi et al. 1999) satellites, shows 
remarkable X-ray variability correlated with strong activ- 
ity at TeV energies (George et al. 1988) on a time-scale of 
~ 10 4 s (Maraschi et al. 1999) and with a flux ratio / — 2. 

X-ray observations of the nearby MKN 766, at z = 
0.013, have been performed by the XMM-Newton satel- 
lite (Boiler, Keill, Tritmper et al. 2000). These observa- 
tions have revealed the presence of a strong X-ray peri- 
odic signal with frequency ~ 2.4 x 10~ 4 Hz and flux ratio 
/~1.3. 

Based on the assumption that the periodic observed 
light curve of MKN 501 is related to the presence of a 
binary system of black holes (one of which emits a jet), 
Rieger and Mannheim (2000) have proposed a method to 
determine the physical parameters of the binary system 
from the observed quantities, i. e. signal periodicity, flux 
ratio between maximum and minimum signal and power 
law spectral index. They adopted the simplifying assump- 
tions that i) the two massive black holes move on circular 
orbits around the common center of mass and ii) the bi- 
nary separation a is such that the gas dynamical time-scale 
T gas is equal to the time-scale for gravitational radiation 
T gw (see also Begelman, Blandford and Rees 1980). 

However, if binary black holes originated from galactic 
mergers, they could be on eccentric orbits and eccentricity 
values up to 0.8-0.9 are not necessarily too extreme (Fitch- 
ett 1987). Of course, due to gravitational wave emission, 
orbits tend to circularize but this happens within a time- 
scale of the same order of magnitude as the merging time- 
scale (Peters 1964, Fitchett 1987). Therefore, if a massive 
binary black hole is found at the center of a galaxy, it 
may happen that the constituting black holes are still on 
eccentric orbits, in which case the method proposed by 
Rieger and Mannheim (2000) does not hold. 

As far as point ii) is concerned, we note that the as- 
sumption that nowadays T gas is equal to T gw is arbitrary. 
Indeed, it is possible that T gas = T gw at some time in 
the past, but since then the binary system evolution was 
driven by gravitational wave emission which shrunk the 
orbital separation faster (see Section 2). 

The aim of the present paper is to generalize the model 
proposed by Rieger and Mannheim (2000) by taking into 
account the effects of the orbital eccentricity of the binary 
system and considering the binary separation a as a free 
model parameter. This is done in Section 2, where the 
model is also applied to three well-observed MKN objects 
(MKN 501, MKN 421 and MKN 766) obtaining a set of 
possible values for the free model parameters jt>, a and e 



corresponding to different values for the masses m and M 
of the two black holes. In Section 3 we study the gravita- 
tional waves emitted by the binary systems at the center 
of the above-quoted MKN objects and analyze the possi- 
bility of detecting them by the LISA gravitational wave 
interferometer. Our conclusions are presented in Section 
4. 

2. Binary black hole model for BL Lacertae 
objects 

In the commonly-accepted evolutionary scenario, giant el- 
liptical galaxies are believed to be the product of mergers 
between spiral galaxies. Since each galaxy likely contains 
a central black hole with mass 10 6 -10 9 M Q (Richstone, 
Ajhar and Bender 1998), merging would sometimes pro- 
duce binary black hole systems separated by a typical dis- 
tance of 0.1 — 1 pc (Begelman, Blandford and Rees 1980)). 

We assume that binary black holes exist in the center 
of many MKN objects and in particular in those of MKN 
501, MKN 421 and MKN 766. As in Rieger and Mannheim 
(2000), the periodicity in the flaring state, observed to- 
wards these objects, is assumed to be the consequence of 
the orbital motion of a relativistic jet in the binary black 
hole. Therefore, the observed signal periodicity has a ge- 
ometrical origin, being a consequence of Doppler-shifted 
modulation. Let us further assume that the jet is emitted 
by the less massive black hole and that the nonthermal 
X — ray /7-ray radiation propagates outwards from the 
core along the jet with Lorentz gamma factor 7?,. The ob- 
served flux modulation due to Doppler boosting can be 
written as 

S(„) = S 3+a S'(u) , (1) 

where a is the source spectral index 1 and the Doppler 
factor is given by 

s y/l-(v*+V? a )/<? 

1 — (v z cos i + V[ s sin 2) fc 

Here, v z is the outflow velocity in the direction of the total 
angular momentum, i is the inclination angle between the 
jet axis and the line of sight and vi s is the component of 
the less massive black hole velocity along the line of sight. 

Let M and m represent the masses of the primary and 
the secondary black hole, respectively. Conventionally, the 
problem of predicting the motion of a system of two grav- 
itationally interacting point masses is simplified by con- 
sidering the equivalent system. This consists of a mass 
M t = M + m fixed in space and acting on a reduced mass 
= Mm/(M + m) that orbits around it. We assume that 

1 Values of the power law index a for the three MKNs of 
interest are found to be 1.2, 1.7 and 2.11 for MKN 501, MKN 
421 and MKN 766, respectively. For more details see Rieger and 
Mannheim (2000), Guainazzi, Vacanti, Maizia et al. (1999) and 
Boiler, keill, Trwmper et al. (2000). 
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the reduced mass moves around the total mass on an el- 
liptic orbit with semi-major axis a and eccentricity e. Let 
(r, 8) be the polar coordinates of /x relative to Mj, so that 
the equation of the orbit is given by 



q(l-e 2 ) 
1 + e cos 8 



(3) 



It is straightforward to demonstrate that the component 
of the reduced mass velocity along the line of sight is given 
by (Smart 1977) 



I G(M + m) 
a(l-e 2 ) 



smt 



(4) 



The component of the less massive black hole velocity 
along the line of sight is thus given by 



M 



vi 



+ m 

which can be rewritten as 



Vis 



Ra^a 



(l-e 2 )i ' 



(5) 



(6) 



where R a = Ma/(M + m) and Sl a = [G(M + m)/a 3 ] 1/2 
is the Keplerian orbital frequency if the mass m moves 
on a circular orbit of radius a around M. The velocity 
vis is maximal for 8 = tt/2 and minimal for 8 = 3tt/2, 
corresponding through equation (2) to the two possible 
values S max and S m i n of the Doppler factor. 

With the assumption that the periodicity in the ob- 
served signal is due to the orbital motion of the binary 
black hole, from equation (1) one obtains the condition 
5max/8min — / 1 /( 3+ "\ where / is the observed maximum 
to minimum flux ratio. Consequently, by using equations 
(2) and (6) we have 



Jl/3+c 



1 ( 1 



/1/3+a + I y s [ n i 



coti (1 - e 2 ) = c , (7) 



which, in the limit e — > 0, reduces to equation (4) in Rieger 
and Mannheim (2000). 

As a consequence of the binary orbital revolution 
around the center of mass, the observed signal period P & s 
is related to the keplerian period P/. = 2it/Q, a by (Rieger 
and Mannheim, 2000) 



Pobs = (1 + z) (l - ^7 cosi) P h , 



(8) 



where i ~ l/^ b and v z ~ c(l - 1/7 6 2 ) 1/2 (Spada 1999). 
From equations (7) and (8), one obtains the following re- 
lation 



M 



dV3 
obs 



(m + M) 2 / 3 ^(l + ^^sim 



/ 1/3+ " ~ 1 (, V Z X V3 

___^l__ cos ,j (i_ e) 



(9) 



2U/2 



which, in the limit e — > 0, reduces to equation (8) of 
Rieger and Mannheim (2000). 

Moreover, from equation (8) one has 



m + M 



2tt(1 + z) (1 - ^cos«) 



obs 



a 
~G 



(10) 



where a represents the semi-major axis of the binary sys- 
tem. 

Therefore, the allowed primary and secondary black 
hole masses m and M are obtained by solving simultane- 
ously equations (9) and (10). In these equations jt, e and 
a have to be considered as free model parameters to be 
determined by the observed data (P i, s , / and a) towards 
the three MKK of interest. In this way, we can estimate 
the masses of the black holes supposed to be at the centers 
of MKN 501, MKN 421 and MKN 766. 

We mention that equation (10) is different from the 
corresponding one in Rieger and Mannheim (2000), where 
the binary separation a is obtained by equating the gas dy- 
namical time-scale (Begelman, Blandford and Rees 1980) 
with the time-scale for gravitational radiation (Peters and 
Mathews 1963). Here, as previously stated, we release this 
assumption and consider a as a free fit parameter. 

For the three MKN objects considered, the secondary 
mass to as a function of the primary one M is shown in 
Fig. 1 for selected values of ^t, = 10, 15 and e = 0.5. For 
each binary system, the intersection between lines with 
the same Lorentz factor determines the masses of the two 
black hole components. 

In Fig. 2 (for MKN 501) we show the primary and 
secondary black hole masses as a function of the orbit 
eccentricity for two values of the Lorentz factor 75 = 10, 30 
and a = 8 x 10 16 cm. In this case, the black hole masses 
are in the range 10 6 — 10 9 M . 

3. Gravitational waves from massive black hole 
binaries in Markarian objects and comparison 
with LISA sensitivity 

In this Section, the massive black hole binaries in the 
MKN objects considered in Section 2 are examined as pos- 
sible sources of gravitational radiation. Due to the orbital 
period of the binary systems considered, it is expected 
that they emit most of the gravitational radiation in the 
low-frequency band (v < 10~ 4 Hz). Gravitational waves in 
this band are only accessible to experiments based on the 
Doppler tracking technique of inter-planetary spacecrafts. 
The Laser Interferometer Space Antenna (LISA), consist- 
ing of a constellation of three drag-free spacecrafts at the 
vertices of an ideal equilater triangle with sides of ~ 5 x 10 6 
km (Hiscock, Larson, Routzahn and Kulick 2000), has an 
optimal sensitivity in the frequency range covering the 
band 10" 5 Hz - 3 x 10" 2 Hz (Vecchio 1999). LISA will be 
able to carry out a deep and extensive census of black hole 
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Fig. 2. The dependence of the primary and secondary mass 
on the eccentricity of the orbit is shown for MKN 501 Lorentz 
76 factors jt = 10 (dotted lines) and jb = 30 (solid lines), 
respectively. Here we have set a = 8 x 10 16 cm. Similar plots 
can be obtained for the MKN421 and MKN766. 



Fig. 1. From the upper to the bottom panel, the required mass 
dependence for binary black hole model is shown in the case of 
MKN 501, MKN 421 and MKN 766, respectively. The thick and 
thin lines represent the conditions expressed in equations (9) 
and (10) with Lorentz factor jt = 10 and jt = 15. The binary 
separation a is set to 5 x 10 16 cm, 2 x 10 14 cm and 1 x 10 13 cm 
for MKN501, MKN421 and MKN766, respectively. We have 
set the orbit eccentricity value equal to 0.5. The intersection 
between lines corresponding to the same Lorentz factor gives 
the masses of the black holes in the binary system considered. 



a circular system with the same semimajor axis. Moreover, 
eccentric binaries will also emit gravitational waves in a 
multitude of harmonics of the orbital frequency ojk associ- 
ated with the semimajor axis a. Thus, an eccentric system 
emits gravitational waves at the frequencies ui n = nuik 
(with n = 1,2,3...), whereas circular binaries emit purely 
in the n = 2 mode. Peters and Mathews (1963), found that 
the orbit averaged power emitted in gravitational waves 
by a binary system with eccentricity e in the nth harmonic 
of the orbital frequency is 



populations in the Universe, providing an accurate demog- 
raphy of these objects and their interactions. In particular 
we expect that binaries of massive black holes could be de- 
tectable by LISA even if the emitted gravitational waves 
are characterized by extremely low frequencies (v < 10~ 4 

Hz). This could be the case for the binaries at the cen- 
ter of the three MKN objects considered under particular 
conditions (e.g. for the same values of a, e and 7&). 

Indeed we expect that in a merging process between 
galaxies, the two colliding black holes at their centers 
would eventually move on orbits with non-negligible ec- 
centricity. This fact should have great importance as far 
as the detection of the emitted gravitational waves is con- 
cerned. In fact, a binary system with an eccentric orbit 
emits significantly more power in gravitational waves than 



dE(n,e) _ 32 G 4 M 2 rn 2 (M + m 



dt 



■g(n,e) , 



(11) 



where the function g(n, e) depends on n and e by the 
following relation 



g(n,e) = — {[J„- 2 (ne) - 2eJ„_i(ne) 



-I — J„{ne) + 2eJ n+1 (ne) - J n+2 (ne) 



+ (1 - e 2 ) [J n - 2 (ne) - 2J n (ne) + J n+2 (ne) 



(12) 



+ 3^2 I J «( ne )] 



■}■ 
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Here J n (with n integer) is the usual Bessel function of 
order n th . The characteristic spectral amplitude at the 
frequency w„ can be evaluated as 

, . , 2 / G dE{n,e) 

h(u n ) ~ J - V - ' 1 , 13 

iru) n D V c A at 

where D is the distance of the gravitational wave source 
from Earth. 

Equations (11)-(13) allow us to evaluate the expected 
spectrum of the gravitational waves emitted by the mas- 
sive black hole binaries assumed to be in the center of the 
considered MKN objects. Assuming 75 = 5, Fig. 3 shows 
the expected gravitational wave amplitude h as a function 
of the frequency cj, for different values of the binary orbit 
eccentricity. In each panel of the Figure we give for com- 
parison the LISA sensitivity curve (the oblique solid line) 
for an integration time of 5 years. As can be noted, in sev- 
eral cases the expected gravitational waves emitted by the 
considered binary systems may be detected by LISA satel- 
lite. Eccentricity values up to 0.8 are not necessarily too 
extreme. In fact, if a black hole binary forms in a merging 
process, a very high value of eccentricity is expected. Of 
course, due to the gravitational wave emission the orbit 
tends to circularize but this should happen within a time- 
scale of the same order of the magnitude of the merging 
time-scale (Peters 1964). 

In Fig. 4, the binary semimajor axis a is shown as 
a function of the orbital eccentricity e and for two dif- 
ferent Lorentz factors >y b for MKN 501, MKN 421 and 
MKN 766, respectively. Each point in the 3D scatter 
plot fits the observations and satisfies eqs. (9) and (10). 
Each point on the plot corresponds to a different distance 
d = h(ui n ) /JiLisAiwn) (vertical axis) between the gravi- 
tational wave emitted spectrum and the LISA sensitivity 
curve for an integration time of 5 years (the minimum 
distance between the curves diminishes as d increases). 
In Fig. 4, only models with infd(w„) > 0.1 are shown. 
Models corresponding to d > 1 in the plot have to be con- 
sidered as observable by the LISA satellite. We notice that 
models corresponding to higher values of the parameter d 
have typically higher masses and lower coalescing times 
with respect to those with lower d values. In any case, 
black hole masses are below 10 9 M . 

4. Discussion 

We have considered some astrophysical implications of the 
existence of massive black hole binary systems at the cen- 
ter of MKN objects. These systems are expected to be 
fairly common in the Universe as a consequence of merg- 
ing processes between galaxies. In fact, following the evo- 
lutionary scenario, giant elliptical galaxies are thought to 
be the result of merging between spirals. Since almost each 
galaxy contains a massive black hole in its center, merging 
process would naturally lead to the formation of massive 
black hole binaries. 



Currently, the only method to probe the existence of 
massive binary black holes is to search for double nuclei 
in the galactic centers. However, as pointed out by Rieger 
and Mannheim (2000) and Yu (2001), a periodic behaviour 
in the observed radio, optical, X-ray or 7-ray light curves 
could be considered as the signature of a massive binary 
black hole. 

In this paper, we focused on the possibility of explain- 
ing the intensity variations observed in some MKN objects 
as the consequence of the motion of binary black holes 
around the common center of mass. We have generalized 
the treatment of Rieger and Mannheim (2000) considering 
in detail the effects of the orbital eccentricity. 

At least three MKN objects (MKN 501, MKN 421 and 
MKN 766) reveal a possible periodic behaviour in the ob- 
served signal. In particular, it has been recently addressed 
whether MKN 501 hosts a massive black hole binary sys- 
tem (Rieger and Mannheim 2000). The same likely holds 
for MKN 421 and MKN 766. The model described in Sec- 
tion 2 allows us to estimate the mass of each black hole in 
the binary system (see Fig. 1) as a function of the Lorentz 
factor 7;,, of the orbital radius a and eccentricity e. 

If massive binary black holes really exist in the cen- 
ters of the MKN objects considered above, they should 
also emit low frequency gravitational waves typically with 
cj < 10~ 4 Hz. In Fig. 3 we have presented the gravitational 
wave spectrum emitted by binaries that fit all the obser- 
vational data (flux ratio /, the signal periodicity P b s and 
the source spectral index a). As one can note, for some val- 
ues of a and e one can obtain curves above the sensitivity 
threshold of the LISA detector (we assume an integration 
time of 5 years) . 

It is obvious that binary systems emitting gravitational 
waves detectable by the LISA satellite have short coalesc- 
ing times. For example, in the case of MKN 501 the av- 
erage coalescing time for the models represented in Fig. 
4 (for 7^ = 5) is ~ 10 3 years and tends to decrease for 
higher values of the parameter d. 

However, for signal detection, some caution is needed 
since the LISA sensitivity curve we use in Figs. 3 and 
4 does not take into account the noise from the galactic 
binaries containing two compact objects. Especially im- 
portant are galactic white dwarf binaries (Hils and Ben- 
der 2000, Nelemans, Yungelson and Portegies-Zwart 2001, 
Nelemans, Portegies-Zwart and Verbunt 1999), primordial 
black holes (Hiscock 1998) and white dwarf binaries (His- 
cock, Larson, Routzahn and Kulick 2000) in the galactic 
halo. These binaries would give rise to a confusion noise 
which lies well above LISA'S sensitivity curve in the fre- 
quency range of interest between 10~ 4 Hz and 10~ 3 Hz 
(see e.g. Fig. 3 in Nelemans, Yungelson and Portegies- 
Zwart 2001). This noise would make even more difficult 
the detection of the considered massive black hole binary 
systems so that, in order to have detectable signals, one 




Fig. 3. The amplitude of the gravitational waves emitted by the binary black holes fitting all the observational data is reported as 
a function of the frequency for the three MKN objects considered. The continuous lines from left to right in each panel correspond 
to orbital eccentricity e = 0.5, 0.6, 0.7, 0.8, respectively (for MKN 501 we added a curve corresponding to e = 0.85). The 
chosen orbital radius a is indicated in each panel. The solid oblique line represents the LISA sensitivity with an integration time 
of 5 years. 



would need more extreme conditions (e.g. higher values of 
the parameter d or, equivalently, shorter coalescing times). 

As a final remark we notice that it is possible that 
there are a lot of hard binary black holes at the center of 
galaxies, MKN objects, Seyfert galaxies, etc., and that we 
can observe directly (from the periodicity of lightcurves at 
high energy) only a small fraction of them. If this is so, de- 
tection of gravitational waves from such systems would be 
a means to detect them. It is worth noting that the num- 
ber of intervening galaxies between Earth and a galaxy at 
redshift z is 

TV ~ 6xl0" 3 ( -^-V ( ^-^) [(l+^) 3/2 -l] ,(14) 

V3k pc ; vo.iM pc -v LV ' iA ' 



where R is the typical galaxy radius and N g the average 
number density of galaxies. Therefore, it is expected that 
at least in some cases gravitational lenses, displaced along 
the line of sight of a gravitational wave source, should dis- 
tort the path of the gravitational waves just as they dis- 
tort the path of the observed light. For gravitational waves 
with frequency v < 10~ 4 Hz, the LISA angular resolution 

is below 10~ 2 steradians (Cutler 1998). Consequently, the 
lensed source is always seen as an amplitude amplification 
(i.e. as a microlensing event) and the observed amplitude, 
at the frequency cj„, is (De Paolis, Ingrosso and Nucita 
2001) 

h max (u n ) ~ y r A^h(u n ) (15) 
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Fig. 4. In the 3D scatter plot, the binary semimajor axis a is shown as a function of the orbital eccentricity e for different values 
of Lorentz factor 7j, for MKN 501, MKN 421 and MKN 766, respectively. Each point in the evidenced regions of the a-e plane 
fits the observations and satisfies eqs. (9) and (10). On the vertical axis, the detection parameter d is also shown (for detalis see 
text). 



where A max is the amplification parameter depending on 
the geometry of the microlensing configuration. This am- 
plification could be in some cases of the order of 10 or 
more, thereby increasing the number of detectable binary 
black holes at galactic centers by future gravitational wave 
detectors like LISA. Even more interesting would be the 
effect of the gravitational wave time-delay due to the dif- 
ferent paths traveled by the deflected signals. This time- 
delay, of about one year as for the lensed Quasars, would 
give rise to a frequency difference due to the binary evolu- 
tion with time. This effect could produce a "beating pat- 
tern" in the wave amplitude which could be searched for 
in future gravitational wave experiments. 
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